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Abstract 

We investigate coarsening and persistence in the voter model by introducing 
the quantity Pn{t), defined as the fraction of voters who changed their opin- 
ion n times up to time t. We show that Pn{t) exhibits scahng behavior that 
strongly depends on the dimension as well as on the initial opinion concentra- 
tions. Exact results are obtained for the average number of opinion changes, 
(n), and the autocorrelation function, A{t) = ^( — ~ t^'^/'^ in arbitrary 
dimension d. These exact results are complemented by a mean-field theory, 
heuristic arguments and numerical simulations. For dimensions d > 2, the 
system does not coarsen, and the opinion changes follow a nearly Poissonian 
distribution, in agreement with mean-field theory. For dimensions d <2, the 
distribution is given by a different scaling form, which is characterized by 
nontrivial scaling exponents. For unequal opinion concentrations, an unusual 
situation occurs where different scaling functions correspond to the majority 
and the minority, as well as for even and odd n. 
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I. INTRODUCTION 



The theory of phase separation, or domain coarsening, has undergone a significant devel- 
opment in the last three decades [|l| . The most important finding is that well-defined ordered 
domains arise and grow with time in such a way that the coarsening process exhibits scaling. 
In other words, at the late stages of the evolution the system is characterized by a single 
length scale L{t) that gives a typical linear size of the domains. It is well established, at least 
for systems with a scalar order parameter, that L{t) ~ with n = 1/2 for non- conserved 
dynamics and n = 1/3 for conserved dynamics. For the Ising spin model, Glauber spin-flip 
dynamics exemplifies the former, while Kawasaki spin-exchange dynamics exemplifies the 
later. 

Several important correlation functions exist. One such function, the autocorrelation 
or, equivalently, the two-time equal-space correlation function, A{t), is defined by A(t) = 
(0(r, O)0(r, if:)), where (f){r,t) is the order parameter. Then, scaling implies: A{t) ~ L^^{t), 
with an exponent A [0. The general two-point correlation function, g{r, t) = (0(0, O)0(r, t)), 
can be expressed through A, namely g{r,t) = L^^G{r / L). Exact results for A are known in 
few cases while the bounds d/2 < \ < d with d the spatial dimension were proposed 

by Fisher and Huse 0. For the 0{m) vector model in the m oo limit, A = (i/2 |^. In this 
study, the value \ = d is obtained for the voter model, defined below. This result indicates 
that both the upper and the lower bounds can be realized. 

It should be noted that in most coarsening processes the dynamics does not exhibit 
a qualitative dependence on the temperature T as long as T < Tc At the critical 

temperature T = T^, the dynamics is generally different, and ordered domains usually do 
not occur. However, the correlation length ^{t) exists and grows with time as ^{t) ~ t^^^, 
where z is the dynamical exponent. The correlation length ^(t) should be considered as the 
analog of the domain size L(t) 0,^], and the the exponent A is replaced by Ac defined by 
Ac{t) ~ C.^^"- In the voter model, temperature is absent but since the dynamics is noiseless, 
the voter model dynamics is zero temperature in nature. However, the "critical" temperature 



2 



is also zero. If one introduces noise by allowing environment-independent opinion changes, 
the system does not coarsen (see, e.g., |^). Thus, we will actually establish \c = d for the 
voter model. A general discussion of the conditions under which the equality \ = \c = d 



holds is given by Majumdar and Huse [10 



In this study, we introduce a family of quantities which provides insight into the "history" 
the the coarsening process. We denote by Pn{t) the fraction of voters who changed their 
opinion exactly n times during the time interval (0,/!:). The first of these quantities, -Po(^), 
is equal to the fraction of persistent voters, i.e. voters who did not change their opinion 
up to time t. This quantity has been introduced independently for two equivalent one- 



dimensional models, the Glauber-Ising model and the single-species annihilation process 
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T4 . Furthermore, the corresponding generalizations to arbitrary dimensions were discussed 
T3[| and [|14|, respectively. Derrida et al. [jl5| established the exact asymptotic decay 



of this quantity, Po{t) ~ t as suggested earlier by numerical simulations. Another 
exact result |]T6| estabhshes Poit) ~ L~^, with [3 = 0.175075. . ., for the ID time-dependent 



Ginsburg-Landau equation at zero temperature. 

For the voter model, several quantities such as the one-time and the two-time correlations 
are exactly solvable in arbitrary dimensions 0,|T^ . These correlation functions allow an exact 



calculation of the average number of opinion changes (n) and other interesting quantities. 
Hence, the voter model is a natural starting point for investigation of Pn{t). Although we 
do not obtain the full distribution, most of its features are illuminated by combining the 
above exact results with heuristic arguments and with the mean-field solution. Generally, 
Pn{t) exhibits a scaling behavior. For d > 2, the scaling function is Poissonian and is peaked 
at n = (n), while for d = 1 the distribution is maximal near the origin. Additionally, for 
unequal opinion concentrations different scaling functions for even and odd opinion changes 
are found. Using random walks techniques, we obtain the even and odd scaling functions in 
the limit of an infinitesimal minority concentration. 

The rest of this paper is organized as follows. In Sec. II, we first solve for Pn{t) on a 
complete graph. We then reexpress some exact relationships for the voter model, in arbitrary 



dimension in terms of the distribution Pn{t). Combining with the mean-field solution, 
these exact relationships allow us to guess the scaling form of Pn(t). This guess suggests a 
usual scaling form in one dimension, and a mean-field like sharply peaked distribution for 
d > 1. These predictions are then compared with numerical data in one, two, and three 
dimensions. In Sec. Ill, we describe exact solution of the mean-field equations for the case 
of initially different concentrations. Then we present exact results for the autocorrelation 
function in arbitrary dimension d, and exact results for the fraction of persistent voters Po{t) 
in one dimension. We proceed by investigating the extreme case of infinitesimal minority 
opinion. In this limit, the model is equivalent to a pair of annihilating random walkers who 
are nearest neighbor at t = 0. Simplifying further the problem to the case of a random 
walker near the absorbing boundary we derive a complete analytical solution. Finally, we 
perform numerical simulations for the case of unequal concentrations and confront the results 
to exact predictions. Finally, we give a brief summary in Sec. IV. 



II. EQUAL CONCENTRATIONS 

In this section we first define the voter model. We restrict attention to the symmetric 
case, i.e., equal opinion densities. We start by analyzing the mean-field theory of the model, 
and then obtain several exact results in arbitrary dimensions. We then present a scaling 
ansatz and check it using numerical simulations. 



A. Mean-Field Theory 



We start by defining the voter model [|T^]. Consider an arbitrary lattice and assume that 
each site is occupied by a "voter" who may have one of two opinions, denoted by + and 
— . Each site keeps its opinion during some time interval, distributed exponentially with 
a characteristic time r, set to unity for convenience, and then assumes the opinion of a 
randomly chosen neighboring site. If a site is surrounded by sites with the same opinion, it 
does not change its opinion. Hence, such dynamics are zero-temperature in nature. Noise 



can be introduced by allowing a voter to change its opinion independently of its neighbors. 
However, a voter system with noisy dynamics does not coarsen, and we restrict ourselves 
to the noiseless voter dynamics. These dynamics are so simple that they naturally arise in 
a variety of situations, see e.g. P,p]7[]. An important link is with the Glauber-Ising model: 



In one dimension, and only in ID, the voter dynamics is identical to the Glauber dynamics. 
This equivalence is not restricted to zero temperature, ID noisy voter dynamics is also 
identical to Glauber dynamics at a positive temperature. 

We now consider the voter model dynamics on a mean-field level, by simply treating all 
sites as neighbors. Such a theory is exact on a complete graph. Moreover, it is expected to 
hold in sufficiently large spatial dimensions. We first consider the symmetric case were the 
opinions concentrations, c+ and c_, are equal, and the interesting case of unequal concen- 
trations will be discussed later. The fraction of sites which have changed their opinions n 
times up to time t, evolves according to 

dPn 

IT = ^"-^ " 

with P_i = to ensure dP^/dt = —Pq. The distribution clearly satisfies the normalization 
condition, J2n Pn = ^, and one can verify that Eq. (|I]) conserves this sum. Solving Eq. (|Tp 
subject to the initial condition Pn(0) = Sno, one finds that the opinion change distribution 
function is Poissonian 

Pnit) = ^e-\ (2) 

In particular, the fraction of persistent voters, i.e. voters who did not change their opinion 
up to time t, decreases exponentially, Po{t) = e~*. The probability that a voter has its initial 
opinion at time t is thus Peven = J2nP2n = (1 + e~^*)/2. Asymptotically, this probability 
exponentially approaches the value 1/2, and therefore voters quickly "forget" their initial 
opinion. 

The distribution is peaked around the average (n) = t, and the width of the distribution, 
a, is given by = (n^) — (n) = t. In the limits, t ^ oo, n — > oo, and {n — t)/\/t finite, 
Pnif) approaches a scaling form 



1 ( n — (n)\ , . 

Pn(t) = -$00 ( , (3) 

where the scahng distribution is Gaussian $oo(a^) = (27r)^^/^ exp(— This infinite- 
dimension scahng solution will be compared below to simulation results in three dimensions. 
To summarize, the quantity incorporates many statistical properties of the system 

such as the probability of maintaining the original opinion, the probability of having the 
original opinion, and the average number of opinion changes. 



B. Exact Results 



We now review several relevant known exact results for the voter model in arbitrary 
dimension d and reexpress them in terms of Pn{t). Both the one- and two-body equal-time 
correlation functions [0,|TB|, are exactly solvable on arbitrary lattice in arbitrary dimension. 



It proves useful to formulate the voter model on the language of Ising spins, i.e., a + opinion 
is identified with -|-1 spin and a — opinion with —1 spin. The state of the lattice is described 
by S* = [^k], the spin-fiip rate Wk(5') = W{Sk —S]s_) reads 

W^k(5) = l--5kE'^k+e„ (4) 

with Zd the coordination number. Here the sum in the right-hand side runs over all Zd 
nearest neighbors. It is convenient to rescale the time variable, t tZd/4:. The probability 
distribution P{S, t) satisfies the master equation 

j^P{S, t)=Y: [w^{S'')PiS^, t) - W^iS)P{S, t)] , (5) 

where the state differs from 5* only at the site k. One can then derive a set of differential 
equations for the spin correlation functions {S^ ... 6*1) = J2s Sk. . . SiP{S, t). The single- and 
two-body correlation functions satisfy discrete Laplace equations 0, 

2^(50 =-^.(5k) + E(^k+e,), (6) 



On a simple (hyper) cubic lattice where Zd = 2d, the general solution for the average 
opinion is given by 

W = e-*^E(^i(0))/k-i(t), (7) 
1 

where Ik{x) is the multi-index function, Ik{x) = Yli<j<dlkj{x), and /„ is the modified 
Bessel function. The autocorrelation Ait) = {So{0)So(t)) is of particular interest since it is 
related to the opinion change distribution via the alternating sum A{t) = J2{—^)"'Pn{t)- The 
autocorrelation is found from Eq. A(t) = e~^'^J2i{So{0)Si{0)) I\(t). In the simplest case 
of completely uncorrelated initial opinions, with equal densities of the opposite opinions, 
(5o(0)^i(0)) = 5io, one finds 

^W = E(-im(t) = [e-*/o(t)]', (8) 

n 

and thus asymptotically, A{t) ~ {2TTt)~~'^^'^ . The diffusive nature of the problem (see e.g. 
Eq. d^)) suggests that the correlation length is given by the diffusion scale, ^{t) ~ 
Therefore, the autocorrelation function scales as ^^'^ for arbitrary d, thus implying that the 
exponent Ac is well defined in all dimensions, and equal to Xc = d as claimed previously. 

The average number of opinion changes, (n) = J2n ^Pn, is simply related to the con- 
centration of "active bonds" (neighbors with different opinions) c_|__ = (1 — (S'iS'i+e))/2: 
d{n)/dt = c_|__. Evaluation of the active bonds density gives the following leading asymp- 



totic behavior 18 



c+_(t) 



l/lnt, d = 2, (9) 
const, d > 2. 

Thus, when d < 2, the density of active bonds vanishes for sufficiently long-time, i.e. coars- 
ening takes place in low dimensions. In contrast, for d > 2, single-opinion domains do not 
arise. This is not very surprising since at the critical point well-ordered domains should not 
necessarily form. Following the above discussion, the average number of opinion changes in 
the limit t oo is obtained by integrating c+_. 



' d<2, 
{n)r^lt/lnt, d = 2, (10) 
J, d>2. 

The above results agree with the mean- field results when d > 2. We therefore expect 
that for d > 2, the distribution function P„(t) approaches the Poissonian distribution of 
Eq. Similarly, the fraction of persistent voters Po{t) should decay exponentially in time 
as well. Interestingly, the exact result for the autocorrelation function indicates a subtle 
failure of the mean-field approach concerning the probability that a voter has its initial 
opinion Peven(^) = I]-P2n(i) = (1 + ^(^)) /2. From Eq. one finds that asymptotically 
-Peven(^) — 1/2 ~ whilc the mean-field approach gives Peven(i) — 1/2 ~ e"^*. Hence 

voters have a stronger than exponential memory, even for d > 2. Despite this erroneous 
prediction, the mean-field theory is successful in predicting most features of the opinion 
change distribution function for d > 2. 

C. Scaling Arguments 

We were unable to find the exact P„(t) distribution for < 2, or higher moments such 
as (n^)- However, combining the above results with scaling arguments proves useful. In one 
dimension, the average number of changes scales as (n) ~ ^/t. We assume that this scale 
characterizes the distribution, or in other words, cr ~ {n) ~ Thus, we arrive at the 
scaling form 

P„W^i^4.(^). (U) 

In general we will use the notation $rf for the d-dimensional scaling function. The nontrivial 
decay of the persistent voter density, Pq ~ t^^^^ [jT5[, implies nontrivial divergence of the 
scaling form ^i{z) ~ z^^^'^, in the small argument limit, z — > 0. The tail of the distribution 
corresponds to a large number of opinion changes by a specific voter and can be estimated by 
an intuitive argument. Such a voter must reside at the boundary between two single opinion 
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regimes, and must change its mind constantly. The probabihty that such a voter changes 
its mind t times (one time per unit time), can be estimated by Pt{t) = exp(— const x t). 
It is also natural to assume that the scaling function rapidly decays for large z, ^i{z) ~ 
exp(— const x z"). Combining this form with Eq. (|TT|) gives Pt{t) ~ exp(— const x and 
consequently, a = 2. To summarize, the limiting behavior of the one- dimensional scaling 
function are 

r z-'/^ 2 < 1; 

~ <^ (12) 
[ exp(— const x z'^) z ^ 1. 

This scaling behavior is different in nature than the scaling behavior for dimensions d > 2. 
While for d >2, a well-defined peak in the distribution function occurs near the average, the 
one dimensional distribution is a decreasing function of n. Moreover, the Gaussian function 
$oo(-2) = (27r)^^/^ exp(— z^) is symmetric around its average, while no such symmetry occurs 
for d = 1, as the distribution is peaked near n = 0. Despite these differences, the tail of 
the distribution of Eq. (p!2| ) agrees with the mean-field distribution of Eq. (^). In fact, the 
above heuristic argument for the large-n behavior is valid in arbitrary dimensions. 



D. Simulation Results 

We implement the voter model using a simple Monte-Carlo simulation. A simple (hy- 
per)cubic lattice is chosen with a linear size L, and periodic boundary conditions are imposed. 
A simulation step consists of choosing randomly an active bond {i.e., a bond between neigh- 
bors with different opinions) and changing the opinion of one of the two voters. After each 
such step, time is incremented by the inverse number of active bonds and the active bond 
list is updated. This simulation procedure is efficient for spatial dimensions d <2 since the 
system coarsens and the number of active bonds decreases as the simulation proceeds. The 
results below correspond to one realization on a lattice of linear dimension L = 10^, 10^, 
and 2 x 10^ in ID, 2D, and 3D, respectively. 

In one dimension, the numerical results confirm the scaling ansatz of Eq. ([TT|), as shown in 
Fig. 1. Interestingly, the maximum of the distribution occurs at n = 1, and the distribution 
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decays monotonously for n > 1. The postulated limiting behaviors of the scaling distribution 
are also confirmed. To test the validity of the mean-field theory, we performed numerical 
simulations in three dimensions. The resulting P„(t) distribution agrees with the Poissonian 
distribution of Eq. (|^), and furthermore the fraction of persistent voters decays exponentially. 
These results indicate that the above mentioned discrepancy regarding the autocorrelation 
function is not crucial in understanding the opinion change distribution function. 

The marginal two-dimensional case is especially interesting. While it is expected that the 
distribution will be roughly Poissonian, some deviation from the mean-field predictions are 
expected. We find numerically that the distribution obeys the scaling form of Eq. and 
exhibits a well defined peak in the vicinity of (n) ~ t/\nt. However, in 2D, the system still 
coarsens, and the distribution exhibits some low- dimensional features. The distribution 
is not a symmetric function of the variable n — (n) (Fig. 2). Additionally, an intriguing 
behavior for the fraction of persistent voters is found numerically (Fig. 3), 

Po{t) ~ exp(-const x In^ t). (13) 

Thus the fraction of persistent voters decays faster than the ID power-law behavior, and 
slower than the mean-field exponential behavior and a naive logarithmic correction to the 
mean-field behavior does not hold. The width of the distribution a is not given by a loga- 
rithmic correction to the mean-field (where a ~ t^^'^), and our best fit gives a ^ t/ In"" t with 
a > 1. The 2D case is difficult from a numerical point of view since logarithmic corrections 
occur, and a large temporal range is required to distinguish such slowly varying corrections 



from algebraic behavior with small exponents |18 . 



III. UNEQUAL CONCENTRATIONS 

Our previous exposition has assumed that the initial concentrations of dissimilar species 
are equal, c+ = c_ = 1/2. The case of unequal concentrations, c_|_ 7^ c_, is interesting as well. 
The reason is that the voter model dynamics has a remarkable feature: Although locally the 
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opinion does change (the dynamics is non-conserved in nature), globally both opinions are 
conserved. On the language of the Ising model it can be said that at zero temperature the 
magnetization remains constant. This hidden integral leads to several peculiarities which 
will be illuminated first on a mean-field level and then for a special case in one dimension. 



A. Mean-Field Theory 

To write the general mean-field theory, it is necessary to distinguish between voters 
according to their initial opinions. Hence, we introduce -P^(t) (P~(t)), the fraction of voters 
with the + (— ) initial opinion that have changed their opinion n times up to time t. If all 
sites are neighbors, these distributions evolve according to 



dt 

dP2n 

dt 



2 (c-P2n-l ~ C^P2n) 



= 2 (c+P^- - c^A-^,) , (14) 

with = 0. The initial conditions are P^(t = 0) = c±Sno- These rate equations reduce 
to Eq. (|1]) for the symmetric case c+ = c_ = 1/2. 

It is again useful to consider the fraction of voters that have (do not have) their initial 
opinion, P^^^ = EnPt (^odd = En ^t+i)- Summation of Eqs. ([T|) gives 

"-^even _ ^f£odd _ 9 2 9 p± /-i c\ 

~ ~ ^""^ ~ ^^^^ 

One can find that the global opinion concentrations c± = P^^n + -^odd conserved and 
that Pq^j = Pq^j. Solving these last rate equations subject to the proper initial conditions 
gives 

^even = C± (c± + C^C-^*) , 

^odd = c+c_ (1 - . (16) 
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The autocorrelation function is then given by 

= (c+ - c_)2 + 4c+c_ (17) 

A voter quickly forgets its initial opinion even if statistically it is more likely to have its 
initial opinion since + > 2c^C-. 

The fraction of persistent voters is found by solving dP^/dt = —2c^P^ and it is found 
that 

p± = c± e-^'^T*. (18) 

Thus, the fraction of persistent voters decays exponentially as well. The decay constant is 
simply given by the density of opposite opinion. This result indicates that even in the case 
of a small concentration of one opinion, the fraction of persistent majority voters decays 
exponentially with time. 

To solve Eqs. (|T^ we introduce the generating functions 



n=0 

oo 



n=0 



This reduces the infinite set of rate equations to four equations 



dF^ 



± 



dt 

^ = 2 (c^^F,±_ - c^F^,,) . (20) 

Expressing F^^ via -F^en; reduce the system of first-order differential equations (^OD to a 
pair of second-order equations for F^^^(t, w) and F^^^{t, w). Solving these equations subject 
to the proper boundary conditions yields 

FLnit, w) = c±e-* fcosh(tA) ± (c+ - ^ (21) 
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where a shorthand notation, A = y 1 — 4c+c_(l — w^), has been used. In principle, one 
then finds Pnit) and Pn{t) by expanding the generating functions. This leads to rather 
cumbersome results. However, the most interesting scaling results correspond to the limit 
t oo, 1 — +0 with (1 — w)t kept finite. In this scaling limit, 1 — A ^ 4c+c_(l — w)t. 
Substituting this into Eq. ( ^Tl) we find -Feven — c^exp(— 4c+c_(l — w)t). Then we find 
F^^, note that in the scaling limit the generating functions become the Laplace transforms 
of Pn{t) and P^{t), and perform the inverse transformation. Finally, we arrive at the 
following scahng results 

Pin ^271 P2n+1 -^2™+! 

2 — 2 — — — 

C_|_ C_ C-|_C_ C-|_C— 

1 

v^27rc+c_t 

In particular, we see that for c+ 7^ c_ the distribution function for the even number of 
changes, P2n = P-m + P2n: is larger than the distribution function for the odd number of 
changes, P2n+i = -P^n+i + An.+i- Eq. ( ^21) suggests that it is possible to avoid these "even- 
odd oscillations", by making a transformation to a modified opinion change distribution 
P„ = Pn + Pn+i- We also note that the scaling distribution in the right hand side of 
Eq. (p^) is identical with the infinite-dimension scaling function, previously obtained for the 
symmetric case. 

B. Exact Results 

Although the above results were obtained using mean-field considerations, similar be- 
havior characterizes the exact solution. By generalizing the solution of Eq. (^, the autocor- 
relation function is found 

Ait) = = (c+ - c^? + 4c+c_[/o(t)e-f . (23) 

The limiting value of the autocorrelation function, (c+ — c_)^, is identical with the mean-field 
theory Eq. (0). Again the conclusion remains the same, at the late stages of the process a 
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(n 



2c+c_ty 



2c4-C-t 



(22) 



single voter opinion cannot be used to determine its initial opinion. Similar to the symmetric 
case, the autocorrelation function decays algebraically rather than exponentially with time. 
Since Peven = (1 + ^(^))/2 > -Podd = (1 " ^(^))/2, wc also learn that a voter is more likely 
to have its initial opinion. 

Mean-field theory suggests that the fraction of persistent voters decays faster for the 
minority. It is interesting to investigate the same for the one-dimensional situation. It 
is instructive to start with the special case of c+ = 1/3 and c_ = 2/3. Let us formally 
split the — opinion into two equivalent sub-opinions. Hence, there are three equiproba- 
ble opinions, one -|- opinion and two — sub-opinions. We now identify this system as the 
zero-temperature three states Potts model, or as a voter model with three opinions. The 
dynamics is unchanged, a voter chooses a nearest neighbor randomly, and assumes its opin- 
ion. Eventually, we will not distinguish between the — sub-opinions. For the kinetic g-state 
Potts model with T = 0, the fraction of persistent spins decays according to Po{t) ~ t~'^^'^\ 
with P{q) = 27r'^ cos^^{^/2q~^ — l/V^) — 1/8, see ||T5[. Indeed, for the symmetric voter 
model, q = 2 and j3{2) = 3/8. The concentration of persistent minority species, P(l{t), 
is equal to the fraction of persistent spins in the g-state Potts model with q = 3. Using 
the notation Po^{t) ~ one has /?+ = (3{3) = 0.5379. Of course, /5_ 7^ since 

changes between — sub-opinions should not be counted. The exponent /5_ can be found 
by allowing a non-integer number of opinions, g = l/c_ = 3/2. This formula is found by 
an analytical continuation to arbitrary q of the relation c = 1/g, which clearly holds in the 
equal-concentration case with an integer q. Therefore, q± = l/c±. For the above example, 
c_ = 2/3, g_ = 3/2 implies /?_ = 0.2349. In general, the concentration of persistent voters 
decays algebraically 
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P(f ~ t-^^ with /3± = — cos-\V2c± - I/V2) - -. (24) 

Following Eq. (p!l|), P„(t) can be written in terms of a simple scaling function in one 
dimension. The z ^ behavior refiects the anomalously large number of persistent voters 
found in the system at long times. On the other hand, Eq. (|2^) implies a difference in nature 
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of the scaling functions for sites of initial + and — opinion, Pn{t) = ^i{n/^/t)/{c±^/t). In 
the limit of large z = nj^^ the tail is dominated by Gaussian fluctuations, while in the 
limit z ^ 0, the anomalous decay of Eq. (|2^ ) determines the behavior. Combining these 
two limits, we have 



In the limit of a vanishing minority opinion concentration, c+ 0, one has /5+ ^ 1, and 

/3_ ^ 2c+/7r ^ 0. 

Both the mean-field results and our numerical simulations, to be described in the follow- 
ing, suggest that distribution of even number of changes dominates over its odd counterpart. 
We expect that the above suggested scaling form holds for the even distribution, or equiva- 
lently, for the modified distribution P„ + P„+i. To summarize, the exact form of the fraction 
of persistent voters combined with scaling considerations suggest that different scaling func- 
tions correspond to the minority and the majority opinions. 



For better understanding of the asymmetric case, it is useful to consider the case of an 
infinitesimal concentration of one opinion, c+ — > 0. We naturally restrict ourself to the 
situation where a single -|- voter is placed in a sea of — opinion. Identifying an interface 
between + and — domains with a random walker, an equivalence to two annihilating random 
walkers who are nearest neighbors at t = 0, is established. The distribution P„(t) is thus 
equal to the fraction of sites visited n times by the two walkers. We further simplify the 
problem by considering the fraction of sites visited by a single random walk with a trap 
as one of its nearest neighbors. Although the two problems are not identical, we expect 
that the results are similar in nature and differ only by numerical prefactors. The reason is 
that the distance between the two random walks itself performs a random walk with in the 
vicinity of a trap. 




exp(— const x z^) z^ 1. 



(25) 



C. Infinitesimal concentrations 
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In the limit of a vanishing opinion concentration, c+ 0, the opinion change density 
Pn(t) is equal to zero. However, if we divide Pnif) by the density of the interfaces, c+c_, and 
then go to the limit c+ 0, we obtain a nontrivial distribution, Ymvcj^^Q Pl^ {t) / c+c^ . This 
distribution gives the total number of links crossed n times by the walker; we will denote it 
by Pn{t). 

As said previously, for the symmetric initial conditions, c+ = c_ = 1/2, the scaling 
behavior of the form P„(t) = t^^^^^i {n/\/t^ is expected. However, for asymmetric initial 
conditions, two different scaling forms, even and odd, should appear. In the present extreme 
case, we expect P2n{t) = t^^^^ ^cvcn (n/y/tj and P2n+i{t) = t^^^'^ ^odd {n/\^^ . We learn 
from Eg . (pSf ) that $evcn = ~ near the origin. Hence, the distribution function 
approaches a time independent form: limt^oo Pn{t) ~ n~^. 

These results can be confirmed by considering the analogy to a single random walk near 
a trap. As the walker will ultimately come to the origin with probability one, every link 
{k — l,k), k > 2 will be crossed an even number of time and so the ultimate distribution 
-P2n+i(oo) = for n > 1 (and Pi(oo) = 1 since the link (0,1) is crossed with ultimate 
probability one by the walker). So, in the extreme case we are considering, the even-odd 
oscillations are obvious and pronounced: The asymptotic even values are positive while the 
odd values are zero. 

In order to compute Pn{oo) for n even, we consider the link {k — l,k). The probability 
that the walker starting at x = 1 reaches for the first time x = k, thus crossing the link 
[k — 1, k), is given by p{k) = 1/k [2^]. Then the ultimate probability that the walker will 
go from site x = A; to site x = k — 1, crossing the link {k — l,k) a second time, is one. The 
probability that the walker starting at x = k — 1 will arrive at a; = before crossing the link 
{k — 1, k) again, is given by 1/k. Therefore, is the contribution of the link {k — 1, k) into 
P2(oo), the average number of links crossed twice by the walker. Thus, we have 

^2(oo) = f:^ = C(2)-l = y-l. (26) 

After having crossed the link {k — 1, k) twice, the walker could cross this link again before 
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reaching the adsorbing barrier at x — 0. Any such crossing from the left happens with prob- 
abihty 1 — 1/k, while the next crossing from the right happens with probability one. Thus, 
we arrive at the remarkably simple formula expressing Pn{oo) through the zeta function 

°° / 1 \" 1 

woo)^i:(i--) 



m=0 



n 



(-ir(C(m + 2)-l) (27) 



For large n, the sum can be approximated by the integral 

1-1/2 1 _ 0-n 1 

P2n{oc) ~ / (1 - 0"-^ = - (28) 

Jo n n 

which confirms the above prediction. 

To determine the scaling functions $cvcn(-2) and $odd(-2), it proves useful to consider 
Pn{x,t), the probability that the walker passes n times through x during the time interval 
(0,i). The Pn{t) is then given by 

^ POO 

Pn{t) = E Pn{x^) ^ / dxPn{x,t). (29) 

x=2 

In this equation and in the following we will treat a; as a continuous variable; in the long-time 
limit, this should be asymptotically correct. 

We then write for Pn{x,t): 

rt rt—ti rt—ti—t2 

P2n{x,t) = / dtipi{x,ti) / dt2P2{h) I dt3P3{x,t3) 

Jo Jo Jo 

E,<2n-i ^^2„p2(^2n) P4 (30) 

and 

rt rt—ti /■t—t\—t2 

P2n+l{x,t) = / dtiPi{x,ti) / dt2P2{t2) / dtsPaix,^) 

Jo Jo Jo 

...j^ dt2n+lP2.{x,t2n)Pb{t- Ylt^ ■ (^l) 

We consider a walker starting at Uq = 1; pi{x,t) is the probability that this walker reaches 
y — X at time t without going to the origin y — 0; P2{t) is the probability that this walker 
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first reaches the origin at time t; P3{x, t) is the probabihty that this walker first passes at the 
origin at time t without passing through y = x; p4^{x,t) is the probabihty that this walker 
with an absorbing boundary aX y = x does not pass through the origin up to time t and 
P5{t) is the probability that this walker does not reach the origin up to time t. Eq. ( pOD 
is cumbersome in form but simple in nature: The formula for P2n{x,t) is just a finite- 
time generalization of Eq. (p7|), namely it corresponds to the situation when a walker has 
performed n oscillations around the link (x — 1, x), and at time t a walker, or his remains, is 
to the left of x. Eq. (^) has been constructed similarly and describes the situation with a 
walker to the right of x at time t. The convolution structure of Eqs. (pOf ) and (|3T| ) suggests 
to apply the Laplace transform. Indeed, Pn{x,s) = dte~'^^Pn{x,t), satisfy 



P2n{x, S) = pi{x, S) (P2(s))" {p^ix, s))" p^ix, s) (32) 

and 

P2n+lix, S) = pi{x, S) (P2(s))'' (psix, s))"p5(s). (33) 

Fortunately, the probabilities pj have been already computed pT[ ]: 

sh(v^) 



Pl[X,S) 



sh{Xy/s) 



P2{s) = e-^ 

_ , , sh((x — l)\/s) , 
P3{x,s)= \ ly ' 34 

Sh{x^/S) 
~ I X 1-P^{X,S) 

Pa[x,s) = - 



s 

1 -P2{s) 



It is in principle possible now to compute various Pn{t)- For example, the contribution to 
Pi{t) from links with k > 2 is 

1 — , , ^ dx 

2 sh{xy/s) 

sh(v^)ln(cth(v^)) (35) 



P,{s) - Pi(l, s) = shiyTs) I 

S J2 

1 - e-v^ 



,3/2 
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(36) 



where the contribution from the first hnk (0,1) is Pi(l,t) = 1 — l/^/Trt, which gives the 
asymptotic value of 

Pi(t) ^ 1 - l/V^t + (t ^ oo). (37) 

y Ant 

We now turn to determination of the scahng functions. In the long-time limit, {t oo), 
corresponding to (s — > 0), Eq. (^) becomes 

g~nv^ ]^ / 1 X n-1 



P2„(x,s)^ ^U-- ' (38) 

s x"^ V xj 

which then implies 

P2n{s)= rfxP2n(a;,s)~ . (39) 

J2 ns 

Performing the inverse Laplace transform one finds 

P2„(t) = -Erfcf4TV (40) 



n VV4t, 

Indeed the anticipated scaling behavior suggested earlier is confirmed with the scaling func 
tion 



$even(^) = ^-'Erfc(^/2) . (41) 



In particular, the limiting forms are 



For the odd distribution, a similar scaling form is expected: 

P2n+l{t) = \KJ^). (43) 



When (s ^ 0), we can use the naive expansion as previously but we should keep the upper 
limit finite, < s~^/^, since the integrand logarithmically diverges on the upper limit: 
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P2n+l{s) - ^ / — 1 



J 2 X \ X, 

Ei{nV^), (44) 



with the exponential integral Ei{y) = duu~^ exp{—u) . Making use of Eq. (^3|) one gets 
another relation for P2n+i{s), 

roo \ / f 



n 



oo 



Vt 



^^^e-^^UlT), (45) 



with q = n^s. 

Thus we obtain the Laplace transform of the function $qjj(T)/-\/T, 

-^\-^^'^,(T)^t:^r'-l.-K (46) 



Performing the inverse Laplace transform, we get 

*-(^' - ^"*Erfc {' \ , ^ exp . (47) 



Performing asymptotic analysis yields 



n « v^, 



;i^\/?exp(-^), n»v^. 

Notice that in the both limiting cases, P2ri+i(^) ^ P2n(t). 

It proves insightful to compute the moments of even and odd distributions, M'P^^^{t) = 
^„>i(2?2)pP2n(i) and M'^^^it) = En>o(2'^+ l)^-P2n+i(^)- Asymptotically, it is easy to com- 
pute even moments 

4PP f p+i\ 

Ml^M=E,t^l\ E, = —^^. (49) 

Eq. ( ^91) is valid only for p > (when (p — > +0), the prefactor Ep diverges). To determine the 
most interesting zero moment, i.e. the total number of links crossed even times M^^^^{t) = 
P^nit), we use the the Laplace transform of Eq. ( P^D to obtain M^^^^{s) ~ ln(l — 
exp(— y^))/s ~ ln(l/s)/2s and eventually, 
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MLnW~(7 + lnt)/2, (50) 

with 7 = 0.577215 the Euler constant. This result is consistent with a direct summation of 
-P2n(oo) = n^^ up to n = i/t. For negative p, even moments are finite, Mp^^^{oo) = 2P({l—p). 

Odd moments behave similarly, M^^^(t) = OpV^"^. A lengthy computation gives the 
prefactor 

O,^ ^1<1,{1-,')' l^^^, (51) 

Eq. (^T]) is valid for all nonnegative p, and in particular the (average) total number of links 
crossed odd times approaches a surprising constant 

Mi.„ = -^rc^/"(^-^'^=ln2. (52) 



Thus although the odd part of the P„ distribution approaches zero as t — oo, the moments 
remain nontrivial. 

D. Simulation Results 

To test the above predictions we performed numerical simulations of the voter model with 
different initial concentrations, in one dimension. The rich behavior predicted by the mean- 
field and the exact results was confirmed by the simulation results. We studied the fraction 
of persistent voters for the case c+ = 1/3, and we found the decay exponents /?+ = 0.54, 
and /3_ = 0.23 for the minority and the majority opinion, respectively. These values are in 



excellent agreement with Eq. (p4D. 



We also confirmed that each of the four functions -P^(t) and -P^+i can be rewritten in a 
scaling with the scaling variable n/\/t. The dominance of the even part of the distribution 
P2n > P2n+i, IS uicely demonstrated by Fig. (4) (one realization of a system of 10^ sites) 
and the asymptotics of the even scaling function Eq. ( pSj) are verified. 

We performed also simulations for the extreme case c+ ^ 0, where one site with initial 
opinion + is in a sea of — opinions. As shown above, this problem is equivalent to the 
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average number of times a link is crossed by two annihilating random walkers. We show 
on Fig. (4) the even and odd scaling functions for 10*^ realizations of this system. The 
asymptotic results Eqs. (^2]) found in the simplified problem of one random walker in the 
presence of an absorbing boundary conditions are verified up to numerical prefactors. In 
particular, the even scaling functions of Fig. (4) is found to behave asymptotically {z 0) 
as $cvcn(-2) — 5/(4z) to be compared with ^evcn{z) 1/z of Eq. (ETI). 



IV. SUMMARY 

We have investigated the voter model, one of the simplest models of non-equilibrium 
statistical mechanics with non-conserved dynamics. We have introduced the set of quantities 
Pn{t), defined as the fraction of voters who changed their opinion n times up to time t. The 
distribution P„(t) was shown to exhibit a scaling behavior that strongly depends on the 
dimension of the system and on the opinion concentrations. For d > 2, the system does not 
coarsen, and the distribution is Poissonian. In one-dimension. We have solved for P„(t) in the 
extreme case when the minority opinion is infinitesimal. The case when the minority phase 
occupies a negligible volume has been studied in the classical work |^ for the conserved 
dynamics and has proven very important in the development of the theory of phase ordering 
kinetics |^. It would be very interesting to generalize the extreme-case solution to arbitrary 
d. 

The quantity Pn{t) reflects the history of the coarsening process. Knowledge of this 
distribution enables insight into interesting quantities such as the fraction of consistent or 
"frozen" sites, the fraction of sites with their original opinion, and the average number of 
changes in a site. This study suggests that P„(t) is a tool for investigations of coarsening 
processes in more realistic models. It possible that a Poissonian Pn{t) generally describes 
systems that do not coarsen, while asymmetric distributions which are pronounced near the 
origin correspond to coarsening systems. 
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FIGURE CAPTIONS 

Fig. 1 Scaling for the symmetric case in one-dimension. The quantity t^/^P„(t) is plotted 
versus n/t^/'^ for different times, t = 10^ 10^, 10^ 

Fig. 2 The distribution function Pn{t) versus n at time t = 2000 in 2D. 

Fig. 3 The fraction of persistent voters in 2D, Po{t) versus In^ t. An average over 300 

samples of linear size L — 10^ for c+ = c_ = 1/2 (solid line) and over 50 samples of linear 

size L = 10^ for c+ = l — c_ = l/4 (dashed line). 

Fig. 4 The even and odd distribution functions for different c+ in ID. i^/^P„(t)/(c_|_c_) is 
plotted versus n/t^^^. Different scaling functions correspond to the even (upper curves) and 
the odd (lower curves) parts of the distribution. The solid hues correspond to the case 
c+ = 1/4 for one sample of linear size L — 10^. The dashed lines correspond to the 
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Figure 4 
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